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1. Introduction 

Cassiopeia A is perhaps the best studied young Galactic supernova Remnant (SNR). It is 
nearby (3.4 kpc) [Reed et al. 1995] and young (~ 325 yr) [Thorstensen, Fesen, & van den Bergh 2001]. 
The wealth of data from ground-based observations in the optical, IR, and radio and from space 
in the optical, x-ray, and y-ray allow us to study its morphology and composition in great detail. 
In [Young et al. 2006] we attempted to narrow down the possibiUties for the progenitor star of Cas 
A using comparisons between computational models and different lines of observational evidence. 
We believe the most probable candidate for the progenitor of Cas A is a 15-25 Mq star in an unequal 
mass binary. This will form the basic assumption of our analysis in this paper. 

This work ties into the larger NuGrid collaboration. We will examine the sensitivity of our 
results to our choice of nuclear rates , determining the dominant rates affecting these yields and 
calculating the error in these yields caused by the uncertainty in the rates. We can then compare 
these errors to current errors in the hydrodynamics, using observed supemovae such as Cas A as 
verification tests in the manner discussed here. 

2. Calculations 

These calculations explore four different progenitor models each for a range of explosion 
scenarios. We use a large set of thermally driven ID explosions with varying delays for a star 
of initial mass 23 Mq and a more restricted range of explosions for a 16 and 23 Mq with the 
hydrogen envelope stripped in a case B binary scenario and a 40 Mq that ends its life as a type 
WC/0 after extensive mass loss. We also examine a 3D explosion of the 23 Mq binary progenitor. 

2.1 ID Explosions 

To model collapse and explosion, we use a I -dimensional Lagrangian code developed by Her- 
ant et al. (1994). This code includes 3-flavor neutrino transport using a flux-limited diffusion 
calculation and a coupled set of equations of state to model the wide range of densities in the 
collapse phase [Herant et al. 1994, Fryer et al. 1999a]. It includes a I4-element nuclear network 
[Benz, Thielemann, & Hills 1989] to follow energy generation. 

2.2 3D Explosions 

Our 3-dimensional simulations use the output of the 1-D explosion (23 M© star, 23m-run5) 
when the shock has reached 10^ cm. We map the structure of this explosion into our 3D Smooth 
Particle Hydrodynamics code SNSPH [Fryer et al. 2006a]. To simulate an asymmetric explosion, 
we modify the velocities within each shell based on angular position. The velocities of particles 
within 30° of the z-axis were increased by a factor of 6 and the remaining particles were decreased 
by a factor of 1.2, roughly conserving the explosion energy. We will refer to these as high velocity 
structures (HVS). At these early times in the explosion, much of the energy remains in thermal 
energy, so the total asymmetry in the explosion is not as extreme as our velocity modifications 
suggest. The large velocity asymmetry results in roughly a factor of two spatial asymmetry between 
the axes. In this calculation, we model the explosion using I million SPH particles. 
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2.3 Nucleosynthesis Post-processing 

Nucleosynthesis post-processing was performed with the Bum code [Young & Fryer 2007], 
using a 524 element network terminating at ^^Tc. Reverse rates are calculated from detailed balance 
and allow a smooth transition to a nuclear statistical equilibrium (NSE) solver at T > lO^'^K. 

3. Explosion Geometry 

At 800 seconds we cannot yet compare the simulation to the detailed observed abundances, 
but we can make statements about the gross morphology. We must account for two features. Most 
of the material is in a circular, flattened structure that is close to the plane of the sky. This is 
well-delineated by the fast moving optical knots, which are largely absent from the center of the 
remnant. In the x-ray, we see that higher atomic number species exist preferentially towards the 
center of the remnant [Willingale et al. 2002]. These features can be reproduced at least qualita- 
tively with a bipolar explosion with the axis oriented close to our line of sight. When the HVS 
expand the bubbles "rupture" after leaving the dense remnant, producing a flattened main remnant 
with nucleosynthesis products of higher A progressively closer to the symmetry axis. 

4. ^^Al Production 

[Ouelette et al. 2007] find that for uniformly expanding SN ejecta, the reverse shock caused 
by the initial interaction of ejecta and a protoplanetary disk develops into a bow shock that deflects 
most of the supernova material around the disk, with injection efficiencies of order of a percent. 
The average ejecta composition is therefore probably not the appropriate quantity to add to the disk 
when considering enrichment by ^^Al . An ejecta knot with the typical characteristics of those in 
the Cassiopeia A SNR would have a mostly or completely unmixed composition, and could deliver 
pure ^^Al -rich material to a disk. Dense optical and infrared knots have densities at least three 
orders of magnitude greater than uniformly expanding ejecta. A dense knot is also the most likely 
part of the ejecta to undergo dust condensation. [Ouelette et al. 2008] find that for dust particles 
with a diameter > 0.01/xm the injection efficiencies approach 100%. For generous assumptions of 
the mass of ejected metals, distance from a supernova, and covering factor of a protoplanetary disk 
we estimate that a disk is likely to encounter only of order one knot that contributes significantly 
to its mass. Since ^^Al or its proxies are unlikely to be detectable in a stellar/planetary system, we 
will consider looking for suitable enriching material in supernova remnants. 

4.1 Detectability of 26a1 Proxies 

Can actually detect the characteristic enhancement or ratios of nucleosynthetic products that 
can reliably indicate the presence of ^^Al in supernova ejecta? Si and S abundances turn out 
to be highly diagnostic (Figure |l]) because of the production channel. The burning raises the Si 
abundance from log(X) ~ — 3 to ~ —1.5, primarily in ^^Si . This is the typical abundance in the 
high ^^Al regions. Temperatures are not high enough to produce significant amounts of S, so S/Si 
drops by a factor of about 10. Oxygen burning produces S efficiently at low entropies ('^O ('^O 
,7)^^S ) and at high T transitioning to QSE, where S is thermodynamically favored over Si. As 
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Figure 1: S/Si (top) and ^^Al (bottom) for model 23m-run5. The ratio drops in the regions of highest ^^Al 
production. Other explosions show a similar pattern. 



a result there is only a narrow range of mass in the star where S/Si drops much below 0.5. This 
coincides with the peak mass fraction of ^^Al . 

We find a similar pattern in the 3D explosion to the ID results (Figure ^ In this case, in 
the spherically symmetric region no material is processed at high enough temperatures to produce 
26ai . The HVS have a large region that consists of material that underwent carbon burning in the 
progenitor star at temperatures of ~ 8 x 10*^K and reached peak shock temperatures of between 1 
and 2xlO^K. This is the ^^Al bubble. The ring consists of material that undergoes Ne burning 
in the progenitor (1 — 1.5 x lO^K) and explosive Ne and/or C burning (2.2 — 2.8 x lO^K). As in 
the ID explosions, S/Si is very low and P production is high. The 3D explosion, with its more 
complicated thermodynamic history, reminds us that ID does not tell the whole story. The yield is 
increased by higher temperatures and higher densities, but ^^Al is rapidly destroyed by succeeding 
reactions, if there is a freezeout after the shock in which density drops rapidly, the high production 
can be achieved without subsequent destruction. This is the reason for the higher ^^Al production 
in the sub-explosive C burning bubble than the explosive ring, the opposite of the pattern found in 
ID. Large amounts of '^F are produced under the same conditions. If this excess flourine is not 
destroyed by subsequent burning it preferentially decays by ^^F (7, a)'^N at temperatures above T~ 
1 X lO'^K. Thus we have materialwith high levels of both N and O, Si, and S. These are potentially 
analogous to the Mixed Emission Knots (MEKs) identified by [Fesen 2001] in Cas A. 

Detection of ^^Al -rich material in core-collapse supernova remnants is a favorable situation. 
The characteristically low S/Si will be undiluted by mixing with ISM material in young remnants, 
so we should be able to predict which ejecta knots have high ^^Al in well-resolved cases like Cas 
A. 0-rich knots with high Mg or Na or moderately enhanced Si are good candidates for a combined 
optical/IR survey. The S/Si ratio of these knots would then provide a strong test of whether they 
would contain ^^Al at high abundance. The Mixed Emission Knots in Cas A may also be excellent 
candidates. These knots are unique due to the simultaneous presence of significant emission from 
N and from O and S. The burning stages in a star that produce O and S destroy N very efficiently. It 
has therefore been assumed that these knots are most likely material from different locations in the 
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Figure 2: Al mass per particle for a 1 O^cm thick slice through the x-z plane. The outer bubble of high Al 
abundance has reached C burning temperatures. The inner ring and bubble reach Ne burning conditions. 

progenitor star that are superimposed on line of sight or somehow mixed during the explosion. If 
these clumps are produced by a nucleosynthetic process they are most likely produced in the same 
conditions favorable for ^^Al production. 
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